Coherent resistance of a disordered ID wire: Expressions for all moments and 

evidence for non-Gaussian distribution 
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We study coherent electron transport in a one-dimensional wire with disorder modeled as a chain 
of randomly positioned scatterers. We derive analytical expressions for all statistical moments of 
the wire resistance p. By means of these expressions we show analytically that the distribution P(f) 
of the variable / = ln(l + p) is not exactly Gaussian even in the limit of weak disorder. In a strict 
mathematical sense, this conclusion is found to hold not only for the distribution tails but also for 
the bulk of the distribution P(f). 
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PACS numbers: 73.23.-b 



I. INTRODUCTION 



It is known that a coherent electron wave in a disor- 
dered one-dimensional (ID) wire of infinite length is 
poncntially localized by an arbitrary weak disorder! 
The resistance p of the ID wire of length L should there- 
fore increase with L exponentially. In fact, the resis- 
tance wildly fluctuates from wire to wire in an ensemble 
of macroscopically identical wires (with disorder in each 
wire being microscopically different) and what increases 
exponentially i^ifihe mean resistance and also the "typi- 
cal" resistanceBEl 

It has also become clear that the resistance p is not 
a self- averaged quantity.^ In fact, the resistance fluctua- 
tions are so huge that (i) the resistance dispersion exceeds 
the mean resistance many orders of magnitude, (ii) the 
higher moments of the resistance exceed the mean resis- 
tance even more drastically, and (iii) the mean resistance 
is much larger than the typical one. These features are 
due to the fact that the moments of p are governed by 
extremely high resistances occurring with an extremely 
low (but nonzero) probability. 

To avoid the absence of self-averaging, the distribution 
P(f) of the variable / Si"^^ + p) was studied instead 
of the distribution P(p).Erljl3 : El In contrast to P(p), distri- 
bution P(f) is well localized around the mean value /. 
It is commonly accepted that for long enough wires the 
bulk of the distribution P(f) is described by the Gauss 
function 
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where A 2 = f 2 — f 2 is the variance, while the tails of the 
distribution P(f) are allowed to be non-universal and 
depend on the model of disorder. In the limit of weak 
disorder it is accepted that A 2 = 2/, i.e., that the dis- 
tribution (Q) obeys the single parameter scaling. The 
two-parameter scaling is accepted to appear for strong 
disorder, where A 2 is not an unambiguous function of 



Interesting to note, the authors of Ref. [if] found 
two-parameter scaling also for weak disorder, namely for 
the Anderson ID disorder at certain conditions. 

In this paper we study coherent transport in a ID wire 
with disorder modeled as a chain of randomly positioned 
scatterers. We derive analytically all statistical moments 
of the wire resistance. By means of these moments we 
prove in the limit of long wires, that the distribution P(f) 
always deviates from the Gauss distribution. The form of 
P(f) for / > / is concluded to be nonuniversal (depen- 
dent on the model of disorder) even in the limit of weak 
disorder. In other words, in realistic wires disorder is 
never weak enough for P(f) to be exactly Gaussian. The 
only approximation of our analysis is the phase random- 
ization hypothesis. We confirm its validity by numerical 
simulations. 

In Sec. II we specify two different model of disordered 
ID wire. As a model I we consider the statistical ensem- 
ble of wires with the same number of scatterers in each 
wire, in the model II we let the number of scatterers to 
fluctuate from wire to wire. In Sec. Ill the moments of 
the wire resistance are derived for both models analyti- 
cally assuming the phase randomization hypothesis. This 
hypothesis is verified in Sect. IV by means of numerical 
simulations. In Sec. V we prove that our expressions 
for the resistance moments are not consistent with the 
Gaussian form of P(f) even in the limit of weak disor- 
der. Discussion is given in Sect. VI. 

II. MODEL OF DISORDERED ID WIRE 

We consider a ID wire with disorder represented by 
random potential 



N 



(2) 



where jS(x — Xi) is the (5-shaped impurity potential of 
strength 7, Xi is the i-th impurity position selected at 



2 



random along the wire, and N is the number of impurities 
in the wire. Since the positions Xj are mutually indepen- 
dent, the distances a — Xi + \ — Xi between the neighboring 
impurities follow the distribution P(a) = JVj exp[— jVja], 
where Nj is the ID density of impurities and Nf 1 is the 
mean distance between the neighboring impurities. 

In the following sections we examine two models. In 
the model I we consider the statistical ensemble of wires 
with iV fixed in each wire to its mean value (N). In 
the model II we fix the wire length L and we let N to 
fluctuate from wire to wire according to the distribution 



Q(N) = (N) N e- {N) /N\. 



(3) 



It is easy to show that this distribution follows from the 
distribution P{a). In both models (N) = LNj. 

The wire resistance p (in units h/2e 2 ) is given by the 
Landauer formulacl 



R{e F ) 

T(s F y 



(4) 



where R and T are the reflection and transmission coeffi- 
cients describing the electron tunneling through disorder 
at the Fermi energy. 

Using eq. (|J) we follow a number of previous localiza- 
tion practitioners. Instead of eq. (0) we could use the two- 
terminal resistance p = 1/T = R/T+l, which involves an 
extra term (unity on the right hand side) representing the 
fundamental resistance of contacts. The resistance (^) 
thus represents the resistance of disorder, directly mea- 
surable only by four-probe techniques. The problem is 
that eq. (^) ignores the effect of measurement probesE3El 
We wish to note that this is not a serious problem in our 
case. First, we examine the regime R/T ^> 1, for which 
the two-terminal resistance p — R/T + 1 coincides with 
eq. (H). Second, with p — R/T + 1 we would arrive at 
the same conclusions as with eq. (|]). Third, in princi- 
ple, one can measure R/T indirectly, by measuring the 
two-terminal resistance and then subtracting unity. 

For disorder (^) both R and T can be obtained by 
solving the tunneling problem 



*L. £- + V(x) 
2m dx 2 



with boundary conditions 



* k (x)=£* k (x) (5) 



tf fe (x - L) = t k e ikx , (6) 



V k (x^0) = e lkx +r k e- 1 

where £ = H 2 k 2 /2m is the electron energy, m is the effec- 
tive mass, and r k and t k are the reflection and transmis- 
sion amplitudes. The coefficients R = \r k \ 2 and T = \t k \ 2 
need to be evaluated at the Fermi wave vector k = k F . 

The reflection coefficient of a single ^-barrier is given 
as #j = fl 2 /(k 2 F + Q 2 ), where Q = m-f/h 2 . We fix 

k F = 7.9 x lOW 1 (7) 

and m = 0.067mo, and we parameterize the 5-barrier by 
Ri. We ignore the fluctuations of Ri as well as the spread 
of the impurity potentials. 



III. RESISTANCE MOMENTS 
A. Model I 

We start with derivation of the mean resistance. As- 
sume that we know the reflection coefficient Rn of a spe- 
cific configuration of N randomly positioned impurities. 
If we add to this configuration an extra impurity at po- 
sition xjv+i, we can express Rn+i through Rn and Ri. 
It is useful to express Rn+i in the forma 



R 



N+l 



Rn + Ri — 2\/R N Rj cos < 



R 



N+l 



{I-Rn^I-Rj) 



(8) 



where 4>n is the phase specified below. Writing eq. (J8|) 
in terms of the wire resistance 



and in terms of 
we get 

Pn+i — Aipjv 4 



Pn 



Ai = 



Rn 



1 - R 



N 
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1 - i?/ 



(9) 



(10) 



Ai-1 



(Xf - l)(pjv + p 2 N ) cos ({>n • 
(11) 

The phase 4>n = 2k F a + (f>o, where a = x^ + i — xn is the 
inter-impurity distance, and </>o is the (a-independent) 
phase due to the reflection by the obstacles. Obviously, 



and 



Pi 



Ri 



Ai-1 



l-Ri 



(12) 



(13) 



Note that p2 depends on <j>i, p% depends on 4>2 and 4>\, 
etc., pn+i thus depends on $jy, 4>n-i% ■ ■ ■ > 4 } 2, and 4>\. 

If we assume that a ^S> 2Ti/k F , then ^jv changes rapidly 
with a and fluctuates at random from sample to sample 
as a fluctuates. The ensemble average of pn+i over J, 
inter- impurity distance xm+i — xn then simplifies too 



2jt 



Pn+i 
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d<t>N Pna 



If we average eq. (11) over </>jv, the term oc cosi 
comes zero. If we then average over <f>N-i, ■ ■ ■ > ' 
we obtain the recursion equation 

Pn+i = Aipjv + -(Ai — 1). 



We solve Eq. (|15|) with initial condition ( |12|) and 
the mean resistance 
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(14) 
f>N be- 

(15) 
obtain 

(16) 
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The higher moments can be obtained in the same way. 
The mth power of eq. ( |TT| ) averaged over </>jy formally 
reads 



„m 

Vn+i ~ 



Ai-1 



If we take into account tha* 

2n 



\/(Ai - l)(p N +p 2 N ) COS (j) N 

n (17) 
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(18) 



and 



2tt 



= 0, 



i / (cos 0) 2 ' 



we easy see that eq. (0) takes the form 

m 



(19) 



(20) 



fc=0 



where coefficients a& (m) are polynomial functions of Ai . 
Averaging each p N over (f>N-i, each p%_ x over <pN-2, 
etc., we finally obtain the recursion relation 



(21) 



A general expression for coefficients ak (m) is given in the 
Appendix where we also derive 



a m (m) = 



Ai — \ X\ — 1 cos i 



(22) 



We can also obtain eq. (123) by comparing the right hand 
sides of Eqs. (|l7| ) and (pOf) for pjv — > oo, where they re- 
duce to [Ai — \JX\ — lcos</>] m p^ and a m (m)p r ^, respec- 
tively. 

We solve E q. (^l|) recursively. Suppose that the N- 
dependence of p 1 ^ can be expressed in the form 



Pn = a m (m)X m + . . .a 1 (m)X 1 + a {m) 



(23) 



For m = 1 eq. ( p3| ) coincides with eq. (|l6|). Therefore, 
Ai in cq. ( p3|) coincides with cq. jl0| ) and ai(l) = 1/2, 
a (l) = —1/2. Once we know Ai, ai(l), and a (l), we 
can solve the problem for m = 2 and determine A2, 02(2), 
ai(2), and ao(2) (see the Appendix |b|). Generally, once 
we determine all A& and all coefficients a n {k) for < n < 
fc < m— 1, we can insert expansion (p3f ) into eq. (|2lf ) and 
compare the iV-independent factors at all A 
gives us linear equations 



k< m - This 



(m)Afe = ^ a i( m ) a k(i) 



(24) 



for all ak(m) with k < m and in addition the identity 
A m = a m (m), i.e., 



Am — 



Xi — <J Aj — 1 cos < 



(25) 



As a last step we calculate the coefficient a m (m) with 
help of the initial condition jl2]). In the Appendix [§| 
this procedure is demonstrated in detail for m — 2. The 
result is 
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where 



A 2 = - (3A? - 1) . 



(26) 



(27) 



Parameters Xk characterize the exponential increase of 
p™ with N . Equation ( p5|) expresses A& analytically for 
arbitrary k, for example, for m — 1 and 2 it reproduces 
relations (0) and (p7|), respectively. We present also 
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(31) 



We do not present explicitly complete expressions for mo- 
ments p^ higher than p 2 N . For further purposes we only 
express the leading term of p 1 ^. We see from ( p5| ) that 
Ai < A 2 < ■ • • < A m . Therefore, for large enough N 



Pn 



■ a < 



™( m ) A^ oc A^. 



(32) 



For completeness, we derive also the mean value of the 
variable /. As in Ref. @, we average over all phases 
the variable Jn = hi (1 + pn) and obtain the recursion 
relation /tv+i = — rn(l— Ri)+/n- We solve this equation 
with the condition po = (i.e, with / = 0) and obtain 



N 



-iVln(l-i?/). 



(33) 



i—k 



No simple analytic expressions exist for higher moments 
7™ For details see Refs. |,|,|lu| 



B. Model II 

In the preceding section the number of impurities, N , 
was kept at the same value for each wire in the wire 
ensemble (model I). In this section we let N to fluctuate 
from wire to wire according to the distribution (|^) while 
keeping for each wire the same wire length L (model II) . 
Thus, to obtain the resistance moments for the model 
II we just need to average over the distribution (||) the 



4 



moments obtained in the preceding section. In particular, 



1.2 



e (A m -l)JV/i = e m{m+l)L/i„ 



(34) 



JV=1 



where we define the mth characteristic length £ m as 

i = Ni (A m - 1) 
^ m m(m + l) ' 

From Eqs. (Fm.M) we obtain the mean resistance 



(35) 



and from Eqs. (p^,|3) the 2nd moment 



^2 = I e 6W«2_I e 2L/? 1 + I. 

F 6 2 3 



(36) 



(37) 



The typical resistance is defined as pt = exp/ — 1. We 
average / (eq. ^) over the distribution ^) and obtain 



p t = expi/£ - 1, 



where 



C 1 = -Nj ln(l - R T ) = Nj In 



Ai + 1 



is the electron localization length. For comparison, 



^ l-Ri 



(38) 



(39) 



(40) 



It is easy to showfij that L/£ m can be expressed as an 
unambiguous function of L/£ a L/£\. This means that 
our models exhibit two-parameter scaling. Only if Rj is 
very small, both lengths converge to the same limit 

£ » & » (Nrfz)- 1 , (41) 

However, £ 7^ £i for any nonzero i?/. 

IV. MICROSCOPIC MODELING 

Our derivation of resistance moments relies on the 
phase randomization hypothesis, i.e., on the averag- 
ing (|l4j). This should be justified in the limit a 3> 2tt /hp, 
that means for 1/Nj ^> 2-Kjkp. Now we test the phase 
randomization hypothesis by microscopic modeling. 

In our microscopic model we select disorder as dis- 
cussed ip_Sect. II, solve Eq. da) by the transfer matrix 
method,E£l and obtain from Eq. (]4|) the resistance of a sin- 
gle wire. We repeat this process for a statistical ensemble 
of wires typically involving 10 6 - 10 9 samples. 

In Fig. |l| we present the distribution P{4>) of the vari- 
able <p, where (f> is the phase entering the right hand side 
of eq. (p"l|). The distribution P(<f>) can be accumulated 
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FIG. 1: Distribution P((j>) for various model parameters. 



either within the ensemble of wires with just two ran- 
domly positioned impurities within each wire, or within 
a single wire into which many impurities are positioned 
one by one. Both procedures give the same results. 

In accord with the phase randomization hypothe- 
sis ([14]), for low impurity density Ni (left panel) we see 
that P{4>) fn const = 1/(2tt). Note that the flat distri- 
bution survives for rather large Rj values. On the other 
hand, when Nj is large, it tends to destroy the flatness 
of P{4>) even for very small values of Ri (right panel). 

Results presented in Fig. [i]are consistent with those in 
Fig. where the mean and typical resistances obtained by 
microscopic modeling are presented for various reflection 
coefficients Rj and various densities iVj. For low Nj our 
microscopic data agree well with our analytical results. 
Note that this is the case also for large Rj. However, 
with increasing Nj the agreement deteriorates. 



V. MOMENTS OF THE RESISTANCE IN THE 
LIMIT OF VERY LONG WIRES. 

In the limit of long wires (N) becomes large and only 
the leading term of the moment becomes important. 
From Eqs. (|32j) and (34) one easy obtains 



A 



(N) 



Model I 



p m oc 



(42) 



e W(A m -i) Model II 
From Eq. |2S it is evident that 



In A ,, 



m > 1 



(43) 



In Fig. U the estimate (^3|) is verified numeric ally . Using 
eq. ( [43|) and (N) = LNi we can obtain from ( |4^ ) 



lnp m (L) 
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(44) 



Now we show that the analytical formulae ( pL4| ) are 
not consistent with the assumption that the distribution 
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FIG. 2: Mean resistance (squares, full lines) and typical resis- 
tance (circles, dashed lines) versus the wire length L. Squares 
and circles are the microscopic model results, full lines and 
dashed lines are graphic representation of the formulas ( pq ) 
and (53), respectively. kF is given in Eq. [?]. Parameters TV/ 
and Ri are varied in such way that the localization length is 
the same (f = 2.7 (im) for each panel. The accuracy of Eqs. 
(pm and (pq) deteriorates with increasing Ni. 



P(f) is Gaussian. To see this clearly, let us average the 
mth power of the resistance 



p™(L) = / df P(f)(ef - I)™ 
Jo 



(45) 



over the Gauss distribution ([j]). The result can easy be 
obtained analytically as 



»(£) = £(-i; 



m — k 



k=l 



"' X exp f ^ + kf). (46) 



In the limit L/£ 3> 1 relation (Eq) reduces to 



p m (L) = exp 



Since A 2 oc L, from (|4^)we have 



+ to/ 



lnp m (L) <x to 2 L. 



(47) 



(48) 



In particular, for weak disorder A 2 = 2f = 2L/£ and the 
leading term in the sum ( ]46| ) reads ex e m ( m+1 * )I '/'' . 

If we compare eq. (f48f ) with our analytical results (44), 
we immediately see that relations (^4|) do not approach 
the dependence cx m?L predicted by relation (|4^). Since 
the higher moments of the resistance are mainly governed 
by the distribution P(f) for / > /, the difference between 
the relations (48) and ( [44] ) is a proof that P(f) deviates 
from the Gauss distribution in the model I as well as in 
the model II. 

It is important to note that these deviations are not 
restricted to the distribution tail / 3> /. It is known 




FIG. 3: m-dependence of In A m for various Ri. 

that the tail of the distribution P(f) is non-universal. 
Tl]) we see that pn+i < (Ai 



From Eq. ([LI]) we see that pn+i < (Ai + \/A 2 — l)pN- 
Therefore, in the model I the value of / never exceeds 
the maximum value / ma x given by 



fn 



In Ar + VAl 3 ! 



In[(Ai + l)/2] 



(49) 



Due to this reason, in the model I the distribution P(f) 
drops to zero for / > / max and some deviations from the 
Gauss distribution ([!]) can be expected to appear already 
for / slightly below / max . The same holds also for the 
model II in which N fluctuates so that the difference 
/max (model II) — / max (model I) is of order of v~N. This 
means that the distribution P(f) drops to zero in both 
models if / is large enough. 

However, this sudden drop to zero is not responsible 
for the nonGaussian behavior represented by eq. (|44|). 
To prove this we now show that p m is governed by the / 
values much smaller that / max . We show that the maxi- 
mum of the function P(f) X e m ^ is positioned at / = f m , 
where f m is much smaller than / max . For the Gaussian 
distribution M) we find 



frn = / + TOA 2 
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/ 



(2m + 1). 



(50) 



This ratio depends neither on L nor on Aj.. Note that the 
ratio / max // does not depend on L but it still depends 
on Ai. In the limit of weak disorder (i?j — > 0) we obtain 
Ai w 1 + 2Ri and /max// ~ 2/V"R7 — > oo. It is thus 
evident that f m <^ / ma x at least in the limit of weak 
disorder. 

From (H) we obtain 



\np m+1 (L) _ m + 1 
In p m (L) m 

for the model I, while for the Gaussian distribution 

lnp m+1 (L) _ m + 2 
In p m (L) m 



(51) 



(52) 
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This proves that P(f) deviates the Gaussian distribution 
already for / from the neighbourhood of fi . As discussed 
above, this region is still far from the distribution tail. 

In the model II this deviation fro m the G aussian shape 
is even more pronounced, because p m {L) increases with 
m much faster than the dependence (M7). This means 
that P(f) decreases for / > / much slower than the 
Gaussian distribution. The slower decrease means that 
the deviation from Gaussian is surely not caused by the 
cutoff at / = / max . 
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VI. DISCUSSION AND CONCLUSIONS 

In conclusion, we have presented two simple models 
of disordered wire which allowed us to express analyti- 
cally all moments of the wire resistance. By means of 
these analytical expressions we have succeeded to prove 
analytically the nonGaussian behavior of the distribution 

P(f)- 

Analytical formulae for the resistance moments were 
obtained assuming the phase randomization hypothesis. 
In Sect. IV we have proven numerically that this hypoth- 
esis is indeed valid for small impurity density TV/. This 
means that for small enough iVj our results are exact. 

In fact, in a strict mathematical sense there is no sin- 
gle parameter scaling in the models I and II, because the 
lengths £ and £i always differ from each other. The dif- 
ference between them is very small in the limit of small 
reflection coefficient, Rj <C 1. Then, numerical experi- 
ment is not able do distinguish between £ and £i and the 
single parameter scaling holds to a good approximation 
for the bulk of the distribution. 

If we accept £_= £i as in eq. (f4l|), then the relations 
(|36|), (|37] ) and ( p8| ) agree with those derived within the 
scaling theory of localization!] Note that the relation ( |4l| ) 
is exact only if the second and higher orders of R[ can be 
neglected. The same condition assures the equivalence 
of eq. (ff8l ) with eq. (Q). Indeed, if we expand A m 
[eq. (p5|)] into powers of Rj and neglect all higher powers 
of i?/, we can interpret the obtained "expansion" as the 
first two terms of the Taylor expansion of the exponential 
function, i.e., 



FIG. 4: The ratio A m / exp [m(m + l)Ri] as a function of Ri 
for m =1 - 5. 



A m = 1 + m(m + l)Ri + O(Rj) 



,[m(m+l)flj] 



(53) 



We show in Fig. |I] that the approximation ( |53| ) is very 
good in the limit of very small Ri and small m. However, 
for any Rj we can find such m that the approximation 
( f53| ) is no longer valid. Therefore, relation (Q) does not 
give the correct Ri dependence for higher moments of 
resistance. This proves that the distribution P(f) is not 
Gaussian even for an infinitesimally small Rj. 

The main difference between the presented results and 
those of the scaling theory of localization is that in our 
model we keep the exact Ri dependence of all A m s while 
in the scaling theory only the linear term in Rj is kept. 
To understand this difference more clearly, let us go back 



to the relation (jTj]). We can approximate Ai as in (|5 
and rewrite ([l5| ) as 



Pn+i = Sp + pn + 25p p N . 



(54) 



Equation (Q) is formally identical with the recursion re- 
lation derived in Refs. ||,f7]. However, in these works it 
is supposed that the increment 8p is proportional to the 
increment SL of the wire length. The terms of higher 
order in So can therefore be neglected and the approxi- 
mation (|53| ) becomes exact. This is not the case in our 
model, where Sp = Rj does not depend on the length 
scale and it is not possible to perform the limit 5L — * 0. 
As this limit plays a crucial role in the derivations of SPS 
m Refs. |and|,' it is understandable that our model does 
not provide us with the Gauss distribution of / predicted 
by these derivations .Ell 
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APPENDIX A 

Coefficients ak(m) in Eqs. (^, ^l|) can be obtained as 
follows. Applying expansion (a+b) m = YliLo (T) a " l ~ l ^ 
and considering Eqs. (Il8|) we obtain from eq. (|l7|) the 
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formula 



I X ■ ' — 



int(m/2) m — 2i i 

E EE 

i=0 j=0 n=0 



m\ / m — 2i 



2t 



A? — 1 cos ( 



>N 



J 



Ai-1 



\ m — 2i— ? rn—i—j-\-n 

x Aj p w . 



(Al) 



To express eq. (Al) in the form (p0|), we choose in the 
triple sum of ( Al) all terms with m — i — j + n = k. We 
write all these terms as a single term ak(m)p > ^ [ , where 

int(m/2) m—2i 

ak(m) — Q(k + — m) 0(m — k — j) 

i=0 j=0 

m\ (m — 2i 



\ 2i j \ j J \k + i + j — m 

,'_> -i ) ^ 1 V yn-2i-J 



Aj — 1 cos ( 



, (A2) 



with 9(x > 0) = 1 and Q(x < 0) = 0. To derive eq. (gg) 
we have also regarded the limits < n < i, which give 
the conditions k + i+ j — m>0 and m — k — j > 0. 

For k = m the funct ion 0(to — k — j) gives the only 
solution j = and eq. reduces to 



int(m/2) 

•(•»>= E (*)>! 



m — 2i 



\\ — 1 COS </>AT 



2; 



Equation QAq) is just binomial expansion of eq. (B2| 



(A3) 



APPENDIX B 



Here we derive the TV-dependence of p 2 N . In accord 
with eq. (23), we assume 



p% = a 2 (2)Af + ai(2)Af +a (2) 



JV 



(Bl) 



where the parameters A2, 02(2), ai(2), and ao(2) have to 
be determined while Ai is known. Also known are the 
coefficients ai(l) = 1/2 and ao(l) = —1/2 [compare Eqs. 
(|) and @ with (§3|) for m = l}. 



Combining Eqs. @, @, and (||) for m=2 we obtain 



P 2 N+1 ^a 2 (2)p 2 N + a 1 (2)pW + a (2), (B2) 

where 02(2) = A2, ai(2) = A 2 — Ai, and ag(2) = (Ai — 
l) 2 /4. Inserting Eqs. fll6| ) and (Bl) into eq. (B2) we get 



a2(2)A 2 A 2 v + ai(2)AiAf + a (2) = 

a 2 {2) [a 2 (2)Af + a 1 (2)Af + a (2)] + 

ar(2) [ ai (l) Af + a (l) ] + a (2). (B3) 

Now we compare the A^-independent factors at A^ = 1, 
Xi , and A^ on both sides of eq. (p|). For A^ we obtain 



(B4) 



ao(2) = a 2 (2) a (2) + ai(2) a (l) + a (2), 



whe re the only unknown parameter is ao(2). Thus, eq. 
(B4) immediately gives 



«o(2) = ~ 

Analogously, for A^ we obtain 

Oi(2) Ax = a 2 {2) 01 (2) + ai(2) oi(l), 

so that 



ai(2) = "2- 
Eventually, for A^ we get 

a 2 (2)A 2 = a 2 (2)a 2 (2), 



(B5) 



(B6) 



(B7) 



(B8) 



which leads to the already known [see eq. (|BJ2j)] identity 
A2(2) = 02(2). In order to calculate 02(2) we have to 
insert the condition p™ = into (Bl). We get 



p2 = = a 2 (2) + ai (2) + ao(2), 



so that 



a 2 (2) 



(B9) 
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